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Abstract 

Assuming neutrinos are Majorana particles having non-zero transition- 
magnetic moments, a global analysis of solar neutrino data, together with 
the new KamLAND data, is presented in the resonant spin flavor pre- 
cession (RSFP) framework. We used Wood-Saxon and Gaussian shaped 
magnetic field profiles throughout the entire Sun. For each magnetic field 
profile, allowed regions from solar data combined with new KamLAND 
data are examined at 95% confidence level (CL). We showed that all al- 
lowed regions are in the large mixing angle (LMA) region and remain the 
same as the fiB value is increased for both magnetic field profiles, contrary 
to the Dirac case studied in our previous work. The electron antineutrino 
flux from the Sun is calculated, and via the results obtained several limits 
are set for fiB. 

1. Introduction 

The magnetic moment solution to solar neutrino problem is alternative pos- 
sible mechanism to MSW [1 — 2] neutrino oscillation solutions. In this solution 
motivated by Okun et al (OVV) [3], strong magnetic field in the Sun changes 
neutrino's helicity which has non-zero transition magnetic moment, as it passes 
through a region with the magnetic field. Shortly after this solution was pro- 
posed, resonance spin flavor precession (RSFP) was proposed by Lim and Mar- 
ciano [4]. In the RSFP scenario, there is a combined effect of matter and 
magnetic field in the Sun on neutrino's spin and flavor precession. Thus, the 
neutrino's spin and flavor flip simultaneously. In other words, left-handed elec- 
tron neutrinos convert to right-handed neutrinos of other types: v,,^ — > I'^j.j^ or 

Vth- t^*' Majorana case v^^ is {y^^^'^ with C being the charge- conjugation 
and called 77^. Matter-enhanced spin-flavor precession of solar neutrinos with 
transition magnetic moments for chlorine and gallium experiments was investi- 
gated in detail by Balantekin et al [5] for Majorana and Dirac case. Raghavan 



et al investigated solar antineutrinos and gave the theoretical framework for the 
solar antineutrino flux calculation [6]. In literature, one can find other studies 
about the solar antineutrino and the experimental limit on its flux [7—10]. In 
[7], observation of the antineutrinos originated from neutrinos in the Sim 
discussed and upper limits on the antineutrino flux are obtained. They used 
the angular distribution of positrons emitted in the reaction of the inverse beta 
decay. Obtaining positrons is a signal for the electron antineutrinos. In that 
study, to flnd the upper limit, statistics of Super Kamiokande (SK) and direc- 
tionality of positrons from inverse beta decay were taken into account. In [8], 
another bound on the antineutrino flux was obtained using SK solar data. The 
antineutrinos contribution to the SK background and angular variations of them 
lead us to set an upper bound on the antineutrino flux. To get a more accurate 
limit, positron angular distribution and antineutrino asimetry was investigated. 
Gando et al [9] gave experimental results for the antineutrino flux. Their limits 
is 0.8% (90%Ci) of the Standart Solar Models' neutrino flux. Miranda et al 
[10] discussed spin flavor precession effect on solar neutrinos assumed to be Ma- 
jorana particle. From the KamLAND [11] constraint on the solar antineutrino 
flux, they put a limit on the Majorana neutrino transition magnetic moment. 

In this work, assuming the fact that all neutrinos are Majorana particles, we 
present a global analysis of solar data combined with the new KamLAND data 
[12] in the RSFP scenario. We obtain allowed regions by using standart least 
squares analysis on the oscillation parameter space, Am^ and tan^ 9, space. 
Our results showed that all allowed regions are in the large mixing angle (LMA) 
region and have the same chi-square value as jiB value is increased for both 
magnetic field profiles. After the combined analysis of solar and KamLAND 
data, we examine the electron antineutrino flux at certain IS.w? and tan'^ 
values. 

This paper is divided as fallows: in section 2 we give general information 
on the evolution equation of Majorana neutrinos having transition magnetic 
moment. In section 3 the relevant magnetic field profiles are given. Detailed 
statistical analysis are given in section 4. In section 5 a theoretical framework 
is given for the calculation of the solar antineutrino fiux. Finally, our results 
and conclusion are presented in section 6. 

2. Evolution equation for Majorana neutrinos in the RSFP frame- 
work 

For two generations, Majorana neutrino flavors are , ,T'e . Altough 
Dirac neutrinos have diagonal and off-diagonal magnetic moments, Majorana 
neutrinos have only off-diagonal (transition) magnetic moments. For the Ma- 
jorana neutrinos that propagate through matter and in a transverse magnetic 
field B, the evolution equation in the case of two generations is 
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where H is given by 
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The matter potentials for a neutral unpolarized medium are given as 

Ve it) = ^ i2N, - 7V„) V, = -^Nr, (3) 

where A^e and Nn are electron and neutron number densities in the Sun, re- 
spectively, and the Gf = 1.16636 x 10~^GeV~^ . In addition in the Sun elec- 
tron and neutron number densities are well approximated by Ng ~ 6Af„ ~ 
2.4 X lO26exp(-r/O.O9i?0)cm-3 [13]. 

One can find Ve resonance from the Hamiltonian through the evolution 

equation. The spin-flavor conversion submatrix for this transition is given by 



Ve IJ.*B 

liB ^ cos 261 -V;, 
Thus Ue — > v., resonance condition is 



(4) 



Aw? 

using the aproximation ~ 6A^„. This value for the electron density required 
for Ve v^ resonance of Majorana neutrinos is greater than the corresponding 
value in the Dirac case. The concequences of this difference will be seen in 
the allowed regions behaviour in the neutrino oscillation parameter space. In 
our analysis, we find results numerically by diagonalizing the Hamiltonian in 
equation (2). The relevant method was discussed in [5]. 

3. Magnetic field profile 

In the literature, various magnetic field profiles have been examined from dif- 
ferent aspects[14 — 18]. Here we considered only two as typical for our analysis 
(see figure 1). 

First, we took the magnetic field profile to be a Wood-Saxon shape of the 
form 

^^"■^ = l + exp[lO(r-ii0)/i?0] 

where Bq is the strength of the magnetic field at the center of the Sun. 

The next magnetic field profile used in the analysis is of Gaussian shape. 

4. Statistical analysis 

For completeness, wc repeat the basics of the statistical methods of our previ- 
ous work. In this case, however, the relevant processes are to be made carefully. 
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Figure 1: Magnetic field profiles: (a) Wood-Saxon sliape; (b) Gaussian shape. 

In the literature, there is a common method often called analysis to find the 
values of the neutrino oscillation parameters Am^, tan^ and to calculate the 
confidence levels of allowed regions and the goodness of a fit [19 — 22]. In our 
analysis, we use "covariancc approach" to find the allowed regions mentioned 
above. By this method, one minimizes the least-squares function 

Afoxp 

ii,i2 

where is the inverse of the covariance matrix of experimental and theoretical 
uncertainties, r['^^'^^ is the event rate calculated in the «th experiment and ij^*'"'^ 
is the theoretical event rate for the ith experiment. The indices indicate the 
solar neutrino experiments: i, ii, Z2 = 1, A^oxp with N^^p = 4. 

Details of the expressions for theoretical event rates for all solar neutrino ex- 
periments, chlorine experiments(IIomestake) [23], gallium experiments (SAGE, 
GALLEX, GNO) [24- 26], Super Kamiokande [27] and SNO [28,29] are given 
in [30]. 

For the global analysis, we need xjcamLAND [12, 31, 32] 

(8) 

We took fluxes and cross sections for event rates and error matrices from 
Bahcall [33]. 

5. Solar electron antineutrino flux 
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If neutrinos arc Majorana type, v,, changes to inside the Sun by SFP. After 
Sun, vacuum oscillation yields v^^v^. This process is given schematically: 

Bsun - Vosc - 

Ve ^ Vii ^ 

To find the electron antineutrino flux on Earth, one need to probability of 
z^e — » Fe transition, Pivf^ Vg): 

P{ve — * '^e) = P{i^e ^ ^/ni SFP) X P(X'^ — > Vg, VacuumOsc.) 

where Ve ^ i^/i transition probability Pive v^;SFP) is calculated numer- 
ically from the equation (2). P{p^ V^,; VacuumOsc.) is the well known 
vacuum oscillation probability given as 



P(77^ ^ X7,; VacuumOsc.) = sin^ 2^sin2(^i?) 1 sin^ 26 

AE 2 

where R is the distance between Sun and Earth. 

Electron antineutrinos is detected observing the positron from the + p — > 
e+ + n reaction. Positron event rate is found from 

N = QoJ dEv j ({E,)R{Ev,E,)<^^SE)(y{E)dE 

where Qq is a normalization constant taking into account the number of atoms 
in the fiducial volume of the detector and its live time exposure. Positron energy 
Ee is 

Ee^ E- 1.293MeV 

and visible energy, Ey = Eg + m^. e{Ee) and R{EY,Ee) are the detection 
efficiency and the energy resolution function of detector, respectively. cf{E) is 
the antineutrino cross section and ^j^^ (E) is the electron antineutrino flux given 

by 

6. Results and conclusions 

In our calculations, we assumed that the magnetic field extends over the 
entire Sun for both magnetic field profiles. All of the calculations have been 
performed for both of the magnetic field profiles. Neutrino spectra are taken 
from the Standard Solar Model of Bahcall and his collaborators [33]. All allowed 
regions were calculated at 95% CL in this work. 

First, we considered only solar neutrino data assuming neutrinos are Ma- 
jorana type and examined allowed regions for all solar neutrino experiments 
( Homestake, Gallium, Super-Kamiokande (SK) and SNO ) using covariance 
approach of statistical analysis. These results are shown in Figure 2 in which 
each column and row are for the same experiment and at the same iJ,B value 
respectively(e.g. in the second row at third column, an allowed region for SK 
experiment at jiB = 2 x 10~^ij,gG is seen). 
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Wc displayed the allowed regions from combined solar neutrino experiments 
in Figure 3. In that figure for both magnetic field profiles, two local best fit 
points are shown for /xB values that are 0,2,5,10 x 10~^/XgG. The best fit 
points are at ( tari^O = 1.63 x IQ-^, Am^ = 6.49 x IQ-^eV^ ) for /iS = 
0, 2, 5 X IQ-VbG and ( tari^O = 1.41 x lO'^, Am^ = 6.49 x IQ-^eV^ ) for jiB 
= 10 X 10~^ij,^G, respectively. 

After this step, we examined the new KamLAND data and displayed allowed 
regions from new binned KamLAND data [12] in Figure 4. 

As a next stage, our global analysis combining solar and KamLAND spec- 
trum analysis were shown in Figure 5 for the same fiB values in Figure 2. 
Regions and the best fit point remained the same for those /xB values, such 
that, the best fit point is at tan^e = 0.26 and at Am^ = 8.66 x lO'^eV^. All 
minimum chi-squares also remained the same for those /iS values. 

It was observed that there are no appreciable differences between the results 
for the two magnetic field profiles we used; we have given only the figures for 
Gaussian case. 

In contrast to Dirac case [34] in the Majorana case we investigated in this 
paper, it seems that there is no appreciable magnetic field effect on the allowed 
regions for Majorana neutrinos in our calculations. 

Finally to put an upper limit on /xiJ for Majorana neutrinos, finally, we 
calculated resulting electron antineutrino flux from the Sun and compared with 
the limit obtained by Super Kamiokande. Our results are given in Ve/^v^ (^B) — 

plane for 8-20 MeV visible energy region and at best fit parameters tan^O = 
0.41, Am^ = 7.1 X 10-^eV2[35] in Figure 6. In that figure, the horizantal 
line shows the upper Ve/^v^i^B) flux limit at SK [9]. From that limit, one 
can put a bound on fiB at these values of /S.rn? and tan^O. From our results, 
fiB < 3.2 X lQ~^iigG and fiB < 3.7 x lO'^/U^G for Wood-Saxon and Gaussian 
shaped magnetic fleld profiles, respectively. 
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Figure 2: The allowed regions for the Majorana case of neutrino parameter space 
for each solar neutrino experiment seperately at fj,B = 0, 2, 5, 10 x 10~ and 
at 95% CL. Each column and row are for the same experiment and at the same 
IJ.B value, respectively(e.g. in the second row at third column, an allowed region 
for SK experiment at = 2 x 10~''/UgG is seen) 
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Figure 3: The allowed regions from combined solar neutrino experiments ( chlo- 
rine, all three gallium, SK and SNO experiments) at the same fxB values and 
CL as in figure 2. Stars indicate the local best-fit points. 
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KamLAND Spectrum 766 Ty Data 
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Figure 4: Allowed region from the new KamLAND spectrum at different confi- 
dence levels. Star indicates the best fit point. 
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Figure 5: Allowed regions from the combined solar plus new KamLAND spec- 
trum analysis at different values {^B = 0, 2, 5, 10 x 10~'^ ij.gG)a.nd at 95% 
CL. Regions and the best fit point remined the same for all fxB values. The star 
indicates the best fit point. 
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Figure 6: Expected Ve/^i,^{^B) for jiB values. The horizantal line shows the 
upper Ve/^u^i^B) flux hmit at SK. 



13 



